Independent pair parton interactions-model of hadron interactions 
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A model of independent pair parton interactions is proposed, according to which, hadron in- 
teractions are represented by a set of independent binary parton coUisions. The final multiplicity 
distribution is described by a convolution of the negative binomial distributions in each of the par- 
tonic collisions. As a result, it is given by a weighted sum of negative binomial distributions with 
parameters multiplied by the number of active pairs. Its shape and moments are considered. Ex- 
perimental data on multiplicity distributions in high energy pp processes are well fitted by these 
distributions. Predictions for the CERN Large Hadron Collider and higher energies are presented. 
The difference between e~^e~ and pp processes is discussed. 



I. INTRODUCTION 

Hadron interactions used to be considered as proceed- 
ing via collisions of their constituent partons. In prepar- 
ton times, their role was played by pions, and the one- 
meson exchange model Q dominated. Pions were treated 
as hadron constituents. Their high energy interaction 
produced a ladder of one-pion t-channel exchanges with 
blobs of low energy pion-pion interactions. This is the 
content of the multiperipheral model. These blobs were 
first interpreted as p mesons and later called fireballs 
0, clusters 0] or clans when higher mass objects 
were considered. Multiperipheral dynamics tells us that 
the number of these blobs is distributed according to the 
Poisson law. It was argued that its convolution with the 
distribution of the number of pions produced in each cen- 
ter can lead to a negative binomial distribution (NBD) 
of created particles first introduced in 0. This supposi- 
tion fits experimental data on multiplicity distributions 
of pp reactions at tens of GeV quite well. However, at 
higher energies this fit by a single NBD becomes unsat- 
isfactory. A shoulder appears at high multiplicities. It is 
quite natural to ascribe it to multiple parton-parton col- 
lisions 0,IE0jE3) which could lead, e.g., to two-, three-, 
and so on, ladder formation [ll|, [ij, uM i and /or to dif- 
ferent (soft, hard) types of interactions They 
become increasingly important as the collision energy is 
increased. Better fits are achieved at the expense of a 
larger number of adjustable parameters. 

This shortcoming can be minimized if one assumes that 
each of the high energy binary parton collisions is inde- 
pendent of others proceeding simultaneously. With this 
supposition, the whole process is described as a set of 
independent pair parton interactions (the IPPI model). 
In fact, we assume democracy in sharing the initial en- 
ergy of colliding hadrons among their constituents. The 
effective multiplicity of particles produced by a pair of 
initial partons does not depend on how many other pairs 
interact or on what these interacting partons are (quarks 



or gluons). While parton energies vary widely at a given 
hadron energy, the mean amount of energy involved in 
each parton-parton collision can be comparable. There- 
fore, we hope that this simplification is valid at very high 
energies for such global characteristics as multiplicity dis- 
tributions which result from an average over the whole 
phase space. If necessary, this supposition can be re- 
laxed by introducing the parton distribution functions 
(PDF) . The inclusive distributions would call for a more 
detailed description. The IPPI model does not imply 
that there are no correlations between particles. They 
are intrinsic in each binary collision and in their convo- 
lution. Surely, further correlations between these inter- 
acting pairs of partons, of both dynamical and kinemati- 
cal origin, can be introduced. Nevertheless, the simplest 
model with minimum parameters and its most general 
characteristics such as multiplicity distributions should 
be treated first. 



II. IPPI MODEL AND MULTIPLICITY 
DISTRIBUTIONS 

For multiplicity distributions, we suppose that a 
parton-parton collision gives rise to a negative binomial 
distribution of its products with parameters independent 
of the colliding partons because the energy is equally 
shared between them. Therefore, the resulting distribu- 
tion of independent pair parton interactions is given by 
a sum of convolutions of processes with different number 
of participating pairs of partons weighted according to 
their probabilities. Thus one can write 

Jmax 

P(n;m, /c) = WjPj{n;m, k) 



X! X! n PNBD{np; m, k). 

i=l (np)P=l 



(1) 
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Here, P(n; m, k) is the probability of creating n particles, 
which depends on the parameters of the NBD distribu- 
tion m and fc, rip is the number of particles produced by 
the pth pair, Wj is the probability for the jth pair to be 
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active and jmax is a number of active pairs. Therefore, 
the following equations are valid 



J2 



= 1. 



(2) 



The symbolical notation X](n ) means the convolution of 
NBD expressions subject to the first equation in jSJl, i.e., 
the sum must be taken only over those parton collisions 
whose multiplicities Up sum up to the total number of 
produced particles n. The NBD shape 



P, 



NBD 



{np;m,k) 



rjup + fc) 
r(np + l)r(fc) 



1 



m 



(3) 



is characterized by two parameters m and fc, correspond- 
ing to the mean multiplicity and the dispersion Di of the 
distribution for a single interaction. 



(4) 



It is our supposition that, weighted by the parton distri- 
bution functions, such an interaction at a fixed parton- 
parton energy leads to a NBD. This is based on the suc- 
cess of low energy fits. In multiperipheral-type models 
with a Poisson distribution of created blobs, this would 
imply a gamma distribution of the decay products of 
these blobs. 

For example, the formula explicitly written for 
three active parton pairs is as follows 



P(n\ m, k) = 



T(n + k) 



r(ni -|- fc)r(n — rii + k) 

(1 + fc/m)" (1 + m/kf r(fc) \^'"'r(n + l) ^ r(fc) (1 + m/fc)'= r(ni + l)r(n - ni + 1) 



W2 



E 



+W3 



Er(ni -f k) y—^ 
Tim + 1) ^ 



r2(fc)(l + m/fc) 



r(n2 -I- fc)r(n - ni — ^2 + fc) 



) r(n2 + l)r(n - 711 - ^2 + 1) 

n2=0 



(5) 



Each of the three terms in this sum represents a negative 
binomial distribution because 

^ r(ni + k)V(n - ni + Ik) _ V{k)V{}k)V{n + \)k) 
r{ni l)r(n - m + 1) ^ r((? + l)A:)r(n + l) 

(6) 

This is a general remarkable property of negative bi- 
nomial distributions: their convolutions result again in 
NBD functions with parameters multiplied by the num- 
ber of convolutions. Thus Eq. can be rewritten as 
follows 

P(n;m, fc) = ^ WjPNBoin; jm, jk). (7) 

This is the main equation of the IPPI model. One gets a 
sum of negative binomial distributions with shifted max- 
ima and larger widths for a larger number of collisions. 
No new adjustable parameters appear in the distribu- 
tion for j pairs of colliding partons. All parameters are 
expressed in terms of the products of parameters for a 
single collision and the number of collisions. Both the 
mean multiplicity and dispersion Z)| for the process with 
J active parton pairs are proportional to j. In the total 
multiplicity distribution, the distributions for collisions 
of j pairs of partons, Pnbd, are just weighted with their 
probabilities Wj, which are determined by collision dy- 
namics and, in principle, can be evaluated if some model 
is adopted (e.g., see [ul 111 El)- 

An increase in the number of interacting pairs of par- 
tons in the IPPI model with energy gives rise to more 



probabilities Wj different from zero. Certainly, all the 
parameters Wj,m,k depend on energy. This dependence 
is implied but not shown explicitly in the above formu- 
las. One can hope that at asymptotically high energies 
the probability for j pairs of independent interactions Wj 
is the product of j probabilities wi for one pair, 

Wj^wi- (8) 

From the normalization condition 

jmax jmax 

£ z., = £ wi = 1, (9) 

one can find wi if jmax , which is determined by the max- 
imum number of parton interactions at a given energy, is 
known. In fact, the value of wi ranges between 1 at low 
energies (for jmax=^) and 0.5 at asymptotics where jmax 
tends to infinity. With energy increase, it approaches the 
second value from above, passing through some thresh- 
olds, and is already quite close to it at the present highest 
energies. 

Thus we are left with only two parameters of the 
model, m and fc, which can be found from fits of ex- 
perimental data. The dependence on the number of col- 
lisions jmax and on the probabilities Wj is determined by 
the behavior of the moments of probabilities 

jmax 

Mr=Yl ^^-J"^' (10) 
J = l 
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as explicitly shown in the Appendix for ranks r < 5. In 
particular, the average multiplicity is given by 

(n) = mMi. (11) 

If one assumes some extrapolation of (n) to higher en- 
ergies, it can be used for prediction of the distributions. 
Let us emphasize that the asymptotic behavior of m is 
directly related to that of the mean multiplicity, 



{nas) 



mas ^ jw{ 



2ma 



(12) 



The value of m is usually quite close to the position of the 
maximum of the distribution. Thus the relation (|12|l tells 
us that in the IPPI model the asymptotic mean multi- 
plicity is about twice larger than the location of its max- 
imum determined mainly by a single parton-parton in- 
teraction. One can expect that the asymptotic relation 
for the probabilities lO becomes valid only at energies 
where four or more pairs are already active. This series 
becomes a polynomial at finite energies. In practice, the 
threshold effects should also be taken into account at fi- 
nite energy. They would somewhat suppress Wj at the 
largest j and, correspondingly, enlarge the role of one- 
and two-pair interactions. 

In 13], the energy dependence of the probabilities 
Wj was estimated according to the multiladder exchange 
model but Poisson distributions were used for each of 
the ladders. The probabilities are given by the following 
normalized expressions 



Pj 




where 



= ln(s/soi^), Zj = 



2C7 



i?2 + a'pij 



(14) 



with numerical parameters obtained from fits of the ex- 
perimental data on total and elastic scattering cross sec- 
tions: 7 = 3.64 GeV-2, i?^ = 3.56 GeV-^, C = 1.5, 
A = ap - 1 = 0.08, a'p = 0.25 GeV-^, so =1 GeV^. 

Below, we will use both possibilities © and H13|) in our 
attempts to describe the experimental data. The prob- 
abilities Wj are different for each (see Tabled). In the 
IPPI model they decrease exponentially with increasing 
number of active partons, while in the ladder model they 
are inversely proportional to this number with additional 
suppression at large j due to the term in parentheses in 
(|13|l . This is the result of the modified eikonal approxi- 
mation. 

We show the values Wj for 3-6 pairs calculated accord- 
ing to Eq. ^ in the left-hand side of Tableland accord- 
ing to Eq. H13|) in its right-hand side. These values jmax 
are chosen because they will be used in the comparison 
with experiment. In particular, we shall choose jmax = 3 



TABLE I: The values of Wj according to Eq. Q (left-hand 
side) and Eq. (right-hand side). 
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Wl 


0.544 


0.519 


0.509 


0.504 


0.562 


0.501 


0.450 


0.410 


W2 


0.295 


0.269 


0.259 


0.254 


0.278 


0.255 


0.236 


0.219 


W3 


0.161 


0.140 


0.131 


0.128 


0.160 


0.153 


0.152 


0.147 


•W4, 





0.072 


0.067 


0.065 





0.091 


0.100 


0.104 


W5 








0.034 


0.033 








0.062 


0.073 


We 











0.016 











0.047 



at 300 and 546 GeV, 4 at 1000 and 1800 GeV, 5 at 14 
TeV, and 6 at 100 TeV (see below). 

One can clearly see the difference between the two ap- 
proaches. The value of wi is always larger than 0.5 in the 
IPPI model while it can become less than 0.5 in the lad- 
der model 0,^5] at high energies. In the ladder model, 
Wj depend explicitly on energy (not only on the jmax 
cutoff). We show their values at 546 and 1800 GeV in 
the right-hand side columns of jmax=^ and 4. Those at 
300 and 1000 GeV are larger for wi by about 1% and 
smaller for by about 3%. When the energy increases, 
the processes with a larger number of active pairs play 
a more important role in the ladder approach compared 
to the IPPI model. Thus, the jmax cutoff is also more 
essential there. 

In principle, one can immediately try a two-parameter 
fit of experimental multiplicity distributions using Eq. 
Q if Wj are known. However, the use of their moments 
is preferred as shown below. 



III. MOMENTS OF MULTIPLICITY 
DISTRIBUTIONS 

The shapes of the multiplicity distributions P{n) usu- 
ally look quite complicated. Often, they are better rep- 
resented by their moments, which also contain complete 
information. The easiest way to define them is to intro- 
duce the generating function 



G(z) = £p(n)(l 



(15) 



In what follows, we will use the so-called unnormalized 
factorial and cumulant ICq moments defined according 
to the formulas 



j^q = ^P(n)n(n-l)...(n-g+l) 



diG{z) 



dz'i 



(P \iiG{z) 



■.=0, (16) 



(17) 



They correspondingly determine the total and genuine 
correlations among the particles produced (for more de- 
tails, see 16]). For q — 1, they define the mean mul- 
tiplicity; the second moment is related to the width of 
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the distribution, etc. The factorial moments are evalu- 
ated from experimental data according to their definition 
(tT?)|l . Both Tq and ICq grow, however, extremely fast with 
their ranks. Therefore, it is more convenient to use 
their ratio Hq — ICq/ J-q, where these dependencies partly 
cancel. This ratio is easy to find from iterative formulas: 



The cumulant moments are written as 



9-1 

1-E 



^r(p + i)r(g-p) 



Hq^,^^^, (18) 



once the factorial moments have been evaluated. Thus, 
both Tq and Hq are determined by experimental data 
according to Eqs. and ifTFI) . 

Recall that these ratios appear quite naturally in QCD 
as the solutions of the equations for the generating func- 
tions of multiplicity distributions 0|. Therefore, their 
use is especially informative because one can compare 
experimental results and model calculations with ana- 
lytical QCD predictions for jets in e+e^ annihilation as 
reviewed in 18]. QCD predicts a very specific behavior 
of the Hq moments as functions of ranks q and energy. 
It has been shown that Hq for jets at the present 

energies (SLC, LEP) should oscillate, and this prediction 
has been confirmed by experimental results 20, 21, 2^ . 
The first minimum is located near q = b dXTP energy. At 
higher energies, this minimum moves to larger values of 
g; the oscillations become less pronounced and disappear 
in the asymptotics where H „ = 1/q^. Moreover, these 
oscillations have been found |18Ll20| even for experimen- 
tally studied pp and AA coUisions. In these cases, they 
can be ascribed to the multicomponent structure of the 
processes. Such a structure is incorporated in the IPPI 
model according to Eq. (Q. We will see if it is enough 
to describe experimental data. 

Let us emphasize that Hq moments are very sensitive 
to minute details of the multiplicity distributions and can 
be used to distinguish between different models and ex- 
perimental data. However, one should be warned that 
the amplitudes of the oscillations strongly depend on the 
multiplicity distribution cutoff due to limited experimen- 
tal statistics (or by another reasoning) if the experiment 
is done at rather low multiplicities. There are no cut- 
offs in analytical expressions for Hq. One can control the 
influence of cutoffs by shifting them appropriately. The 
qualitative features persist nevertheless. In what follows, 
we consider very high energy processes where the cutoff 
due to experimental statistics is practically insignificant. 

The IPPI model predicts new special features of the 
moments J-'q and Hg. The factorial moments of the dis- 
tribution ((TJ are 



^9=E 



T{jk + q) ( m\<i 



(m\i „ , (m\i , , 



with 



/9(fc) = E 



r(jfc + q) 



.=1 ^ ^^^^ j=i 

(20) 



^<«(t)'^ 



(21) 



The explicit fc dependence of jq(k) and Kq{k) for g < 5 
is shown in the Appendix. For Hq moments one gets 



H„ 



9-1 

1-E 



^r(p + i)r(g-p) 



Hq^p 



fpfq- 

~17 



(22) 



Note that according to Eq. (|22|) Hq are functions of the 
parameter k only and do not depend on m in the IPPI 
model, because the m dependence of factorial and cumu- 
lant moments is the same. This remarkable property of 
the Hq moments provides an opportunity to fit the same 
results with a smaller number of parameters. If the Wj^s 
are given by Eq. 0, the only adjustable parameter left 
is k. These moments decrease with increase of k and q. 

Once the parameter k is found from fits of Hq, it is 
possible to get another parameter m by rewriting Eq. 
(|T^ as follows: 



•^9 

fq{k) 



1/9 



(23) 



This formula is a sensitive test for the whole approach 
because it states that the definite ratio of g-dependent 
functions to the power 1/q becomes q independent if the 
model is correct. Moreover, this statement should be 
valid only for those values of k that are determined from 
Hq fits. Therefore, it can be considered as a criterion for 
a proper choice of k and for the model validity, in general. 

One substitutes the experimentally determined values 
of factorial moments, divides them by the theoretical 
functions fq given by Eq. H20|) . and examines whether 
this ratio to the power l/q is independent of g at fc val- 
ues found previously from Hq fits. If the answer is posi- 
tive, the parameter m is known according to Eq. (|23l) . If 
not, the model should be modified. With parameters k 
and m found, one can try to fit the shapes of experimen- 
tally measured multiplicity distributions directly. This is 
another test of the self-consistency of the IPPI model. 

At the same time, the value of m determines the po- 
sition of the peak of the multiplicity distribution. For a 
given (n) , it can be used to check the choice of probabil- 
ities Wj according to Eq. ifTTl) . 

Recall that both parameters m and k depend on the 
energy of the colliding hadrons s. This dependence can 
be determined from fits of experimentally found values 
of Hq and J-q as explained above. To extrapolate it to 
higher energies, one should use some guesses. Since m 
has the meaning of the average multiplicity of a binary 
parton collision, it should behave similarly to the mean 
multiplicity of the whole process. The latter is usually 
fitted by a logarithmic dependence with some log-squared 
terms added. No experience has been gained yet for the 
parameter k. 
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The Poisson distribution possesses the same property 
of convolutions which made it possible to get Eq. iQ for 
NBD distributions. Therefore, all the above relations are 
valid for a model with convoluted Poisson distributions. 
Actually, they can be obtained in the limit fc ^ oo. For 
example, the factorial moments are _^^^°*^*°") — m'^Mq. 

IV. COMPARISON WITH EXPERIMENT 

We have compared IPPI model conclusions with ex- 
perimental multiplicity distributions of the E735 Collab- 
oration "23] for pp collisions at energies 300, 546, 1000, 
and 1800 GeV extrapolated 0, [3| to the full phase 
space. The multiplicity of charged particles was divided 
by 2 to get the multiplicity of particles with the same 
charge. Then the above formulas for the moments were 
used. Correspondingly, the parameters m and k refer to 
these distributions. 

An analysis of experimental data done in has 
shown that two parton pairs are already active at en- 
ergies above 120 GeV. The thresholds for triple or more 
parton-parton collisions are less definite. They depend 
on the form of the multiplicity distribution adopted for 
a single collision. We assume that three parton pairs are 
active at 300 and 546 GeV and four at 1000 and 1800 
GeV with NBDs for a single collision. We use these val- 
ues in our calculations. 

Factorial and Hq moments were obtained from exper- 
imental data on P{n) according to Eqs. (|16l) and (|18|l . 
Experimental Hq moments were fitted by Eq. ^ to get 
the parameters k{E) of the IPPI model. We show in Fig. 

1 how perfect are these fits at 1.8 TeV for k equal to 3.7 
(solid line) and 4.4 (dash-dotted line). At this energy, we 
consider four active parton pairs with Wj given by Eq. 
© (the second column in Table HJ. It is surprising that 
the oscillations of Hq moments are so well reproduced 
with one adjustable parameter k. The general tendency 
of this quite complicated oscillatory dependence is clearly 
seen. 

With these values of the parameter fc, we have checked 
whether m is constant as a function of q. Experimental 
factorial moments and IPPI values for fq were inserted 
in Eq. (|23|l . The m{q) dependence is shown in Fig. 2 for 
the same values of k=AA (squares) and 3.7 (circles) and 
for a much larger value 7.5 (triangles). The constancy 
of m is satisfied with an accuracy better than 1.5% for 
k — AA up to g = 16. The upper and lower lines in Fig. 

2 demonstrate clearly that this condition substantially 
bounds the admissible variations of k. 

It is well known that experimental cutoffs of multiplic- 
ity distributions due to the limited statistics of an exper- 
iment can influence the behavior of Hq moments. Conse- 
quently, they impose some limits on the q values that can 
be considered when a comparison is done. Higher rank 
moments can be evaluated if larger multiplicities have 
been measured. To estimate the admissible range of q, 
we use the results obtained in QCD. Characteristic mul- 
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2 4 6 8 10 12 14 16 

q 

FIG. 1: A comparison of Hq moments derived from experi- 
mental data at 1.8 TeV (squares) with their values calculated 
with parameter k — 4.4 (dash-dotted line) and 3.7 (solid line). 

tiplicities that determine the moment of the rank q can 
be found. By inverting this relation, one can write the 
asymptotic expression for the characteristic range of q 
[25j. This provides the bound qmax ~ Crimax/ {n) where 
C ~ 2.5527. However, it underestimates the factorial 
moments. Moreover, the first moment is not properly 
normalized (it becomes equal to 2/C instead of 1). The 
strongly overestimated values (however, with a correct 
normalization of the first moment) are obtained if C is 
replaced by 2. Hence, one can say that the limiting values 
of q are given by the inequalities 

2nrnax/{n) < q,nax < Cn„iax/{n)- (24) 

The ratio n„iax/ (n) measured by the E735 Collaboration 
at 1.8 TeV is about 5. Thus, qmax should be in the 
interval between 10 and 13. The approximate constancy 
of m and proper fits of Hq demonstrated above persist 
to even higher ranks. 

The same-charge multiplicity distribution at 1.8 TeV 
has been fitted with the parameters m = 12.94 and k — 
4.4 as shown in Fig. 3 (solid line). To estimate the accu- 
racy of the fit, we calculated X]n=i[^t("-) ~ Pe{n)]'^/A'^ 
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FIG. 2: The g dependence of m for fc = 4.4 (squares), 3.7 
(circles), and 7.5 (triangles). 



over all 125 experimental points. Here, Pt,Pe are the 
theoretical and experimental distributions and A is the 
total experimental error. It includes both statistical and 
systematical errors. Note that the latter are large at low 
multiplicities in the E735 data. This sum is equal to 50 
for 125 degrees of freedom. No minimization of it was at- 
tempted. This is twice better than the three-parameter 
fit by a generalized NBD considered in (2^]. A Poisson 
distribution of particles in binary collisions is completely 
excluded. This is shown in Fig. 3 by the dash-dotted 
line. 

We show in Fig. 4 the decomposition of the fit in Fig. 
3 to processes with different numbers j of parton pairs 
involved in collision. It is seen that the locations of their 
maxima are approximately proportional to j . 

The same procedure has been applied to data at ener- 
gies 300, 546, and 1000 GeV. As stated above, we have 
assumed that three binary parton collisions are active at 
300 and 546 GeV and four at 1000 GeV. We plot in Figs. 
5 and 6 the energy dependence of the parameters m and 
k. The parameter m increases logarithmically with en- 
ergy. This is expected [see (|11|) ] because increase of Mi 
due to increasing numbers of active pairs at these ener- 



0.06 




100 



FIG. 3: The muhiplicity distribution at 1.8 TeV and its 
fit at m = 12.94, fc = 4.4 (solid line). The dash-dotted line 
demonstrates what would happen if the NBD were replaced 
by Poisson distribution. 



gies leads to a somewhat faster than logarithmic increase 
of the average multiplicity in accordance with experimen- 
tal observations. The energy dependence of k (crosses) 
is more complicated and rather irregular. 

We tried to ascribe the latter to the fact that the ef- 
fective values of fc, which we actually find from these fits, 
depend on the effective number of parton interactions, 
i.e., on the Wj variation at a threshold. The threshold 
effects can be important in this energy region. Then the 
simple relation ((SJ is invalid. This infiuences the func- 
tions fq{k) [Eq. 1(201] and, consequently, Hq calculated 
from Eq. H22|l . One can reduce the effective number of 
active pairs to about 2.5 at 300 GeV and 3.5 at 1000 GeV 
if one chooses the following values of Wj: 0.59, 0.34, and 
0.07 at 300 GeV and 0.54, 0.29, 0.14, and 0.03 at 1000 
TeV instead of those calculated according to Eq. © and 
shown in Table U This gives rise to values of k which are 
not drastically different from the previous ones. How- 
ever, the quality of the fits becomes worse. Fits with two 
active pairs at 300 GeV and three pairs at 1000 GeV fail 
completely. 
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100 1000 10000 100000 

E, GeV 

FIG. 5: The energy dependence of m (squares) and its linear 
extrapolation (circles at 14 and 100 TeV). 



FIG. 6: The values of k as calculated with Wj satisfying the 
relation Q (squares). 



Hence, we have to conclude that this effect results from 
some dynamics of the hadron interactions that is not un- 
derstood yet and should be incorporated in the model. 
The preliminary explanation of this effect could be that 
at the threshold of new pair formation the previous ac- 
tive pairs produce more squeezed multiplicity distribu- 
tions due to the smaller phase-space room available for 
them because of the newcomer. Therefore, the single 
pair dispersion decreases and the k values increase. This 
would imply that thresholds are marked not only by the 
change of Wj shown in Tabled but also by the variation 
of the parameter k. 

The threshold effects become less important at higher 
energies. We assume that there are five active pairs at 14 
TeV and six at 100 TeV. Then we extrapolate to these 
energies. The parameter m becomes equal to 19.2 at 14 
TeV and 25.2 at 100 TeV if logarithmic dependence is 
adopted as shown in Fig. 5 by the straight line. The 
predicted multiplicity distributions are plotted in Fig. 
7. We choose two values of k equal to 4.4 (solid line) 
and 8 (dash-dotted line) for 14 TeV. Low multiplicities 
are suppressed at larger k, and the maximum is slightly 
shifted to higher multiplicities. The shape of the tail is 
practically unchanged. For 100 TeV, we show only the 
prediction for k — 4.4 (dashed line) because increase of 
k leads to the same qualitative effect as for 14 TeV. The 
oscillations of Hq still persist at these energies (see Fig. 
8). The minima are, however, shifted to g = 6 at 14 TeV 
and 7 at 100 TeV as expected. 

The fit at 1.8 TeV with an approximation oiwj accord- 
ing to the ladder model H13|) with a NBD for a binary par- 
ton collision is almost as successful as the fit with values 
of Wj given by the IPPI model. However, some difference 
at 14 TeV between these models is predicted (compare 
the solid and dotted lines in Fig. 7, both obtained for 
k = 4.4). This difference becomes more pronounced at 
100 TeV. To keep the same mean multiplicity in both 
models at the same energy, we have chosen different val- 
ues of m as dictated by Eq. (|ll|l and the Wj values shown 
in TableHl namely, their ratios are mippi /miad =0.988, 
1.039, 1.123, 1.228 iov jmax='i, 4, 5, 6, respectively. This 
shows that the maximum of the distribution moves to 
smaller multiplicities and its width becomes larger in the 
ladder model compared to the IPPI model with energy 
increase. 

Certainly, one should not overestimate the success of 
the IPPI model in its present initial state. It has been ap- 
plied just to multiplicity distributions. For more detailed 
properties, say, rapidity distributions, one would need a 
model for the corresponding features of the one-pair pro- 
cess. Moreover, the screening effect (often described by 
the triple Pomeron vertex) will probably become more 
important at higher energies. All these features are im- 
plemented in some way in the well known Monte Carlo 
programs PYTHIA HERWIG [H and DPM-QGSM 
|Tll 12j. However, for the last one, the multiplicity distri- 
bution for a single ladder is given by the Poisson distri- 
bution of emission centers (resonances) convoluted with 
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FIG. 7: The same-charge muhiphcity distributions at 14 TeV 
and 100 TeV obtained by extrapolation of parameters m and 
k with five active pairs at 14 TeV and six at 100 TeV (for 
the IPPI model: solid line, 14 TeV, k = 4.4; dash-dotted line 
14 TeV, fe = 8; dashed line 100 TeV, k = 4.4; for the ladder 
model: dotted line 14 TeV, k = 4.4). 



their decay properties, and the probabihties Wj contain 
several adjustable parameters. It differs from the IPPI 
model. The present approach proposes more economic 
way with a smaller number of such parameters. Concern- 
ing the further development of event generator codes, it 
is tempting to incorporate there the above approach with 
a negative binomial distribution of particles created by a 
single parton pair, and confront the results with a wider 
set of experimental data. This has not been done yet 
for the IPPI model, and we intend to work on it later to 
learn how it influences other characteristics. 

It would also be interesting to see whether this model 
is valid for AA collisions as well whether the collective 
effects (saturation?) prevent its application there. This 
work is in progress now. 



0.35 




-0.05 - 



-0.1 I 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 

q 

FIG. 8: The behavior of Hq predicted at 14 TeV (fe = 4.4, 
solid line; fe = 8, dash-dotted line). 

V. ARE e+e" AND pp SIMILAR? 

This question was raised by a recent statement of the 
PHOBOS Collaboration that the energy behavior 
of mean multiplicities in all processes is similar. It was 
concluded that the dynamics of all hadronic processes is 
the same. In addition to our general belief in QCD, we 
cannot claim that other characteristics of multiple pro- 
duction processes initiated by different partners coincide. 

At first sight, QCD fits of multiplicity distributions in 
e+e~ collisions and IPPI model fits of pp collisions are 
completely unrelated and cannot be compared. There is, 
however, one definite QCD prediction that allows us to 
ask the question whether QCD and the IPPI model are 
compatible. This is the asymptotic behavior of Hq mo- 
ments in QCD. They should behave [l3 as H^" = 1/q'^. 
One can also determine the asymptotics of Hq moments 
in the IPPI model and compare both approaches. The 
asymptotical values of the probabilities Wj (|HJl and their 
moments Air (|10|l for r < 5 are as follows: 

= 0.5^ Ml = 2, M2 = 6, Ms = 26, 

M4 = 150, Ms = 1082. (25) 
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Inserting them in the expressions for Kg and fq given in 
the Appendix, one can evaluate the asymptotic behavior 
of the Hq moments in the IPPI model at any parame- 
ter k. All asymptotic Hq are decreasing functions of k. 
Their minimum values are reached at k — > oo, i.e., for 
a convolution of Poisson distributions. They are given 
by the ratio of the coefficients in front of the leading fc'? 
terms in Kq and fq (see the Appendix) and are equal to 

tt(p) _ 1 tt(p) _ 3 (P) _ 13 (P) _ 75 

H, -3, H, --, H, --, H, -—. (26) 

These values are noticeably larger than QCD predictions 
of Since they are even larger for any finite param- 

eter k, we have to state that QCD and the IPPI model 
have different asymptotics. In other words, this implies 
that Eq. |j23), considered as an equation for fq with 
Hq = inserted in it, does not have a solution with 
asymptotical values of Mr [Eq. (|25|l ] in the IPPI model. 

It is an open question whether other asymptotic rela- 
tions for Wj different from Eq. JSJ can be found which 
would lead to the same behavior of Hq moments in pp 
and e+e" collisions, i.e., if a solution of Eq. 1)22(1 can be 
found for some values of Mr different from those given by 
Eq. (|25|l . Only then one can hope to declare an analogy 
between these processes. 

Moreover, it has been found from experimental data 
[isl 1^ that the amplitudes of oscillations of Hq mo- 
ments increase for more composite collidin g p articles. 
The anomalous fractal dimensions also differ [1^ becom- 
ing smaller for AA compared with pp and even more with 
e~^e~ . Thus, there is no direct similarity of e~^e~ and pp- 
collisions. 



VI. CONCLUSIONS 

To conclude, a model of independent pair parton inter- 
actions has been proposed. It is assumed that hadronic 



interactions proceed through independent parton-parton 
collisions and each of the binary collisions gives rise to a 
negative binomial distribution of secondary particles with 
the same parameters m and k. The resulting distribution 
is described by a weighted sum of NBDs whose parame- 
ters are equal to the single collision values m and k mul- 
tiplied by the number of pairs. Thus no new adjustable 
parameters appear. Multiple binary parton collisions are 
assumed to become more important as energy increases. 
A comparison with experimental data at 300, 546, 1000, 
and 1800 GeV has shown good agreement. Predictions 
for the CERN Large Hadron Collider and higher ener- 
gies are presented. It is demonstrated that asymptotic 
QCD predictions for e^e^ multiplicity distributions dif- 
fer from the asymptotic results of the IPPI model for pp 
processes. Further work on Monte Carlo implementation 
of this model is in progress. 
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APPENDIX A 



The functions fq{k) and Kq{k) are qth order polynomi- 
als of k with coefficients determined by the moments Mr 
with r < q. Their expressions for q < 5 are as follows: 



fi = ni^ Mik, f2 = M2k^ + Mik, K2 = {M2 - Ml)k^ + Mik, 
/g = Mg/fc^ M'hk^ + 2Mifc, K3 = (A/a - 3M1M2 + 2Ml)k^ + 3(M2 - Aff )fc2 + 2Mik, 

fi = Mik^ + QM^k^ + llM2k^ + 6Mik, 
K4 = (A//4 - 4A'/iAf3 + 12A/f Af2 - 3Af| - 6M^)k* + 6(M3 - 3AfiAf2 + 2M^)k^ + ll(Af2 - Mf)k^ + QM^k, 

/g = Afsfc^ + lOMik^ + SdMsk^ + 5OM2P + 24:Mik, 
K5 = (Afs - 5AfiA/4 + 20Af fMg - 60MfM2 + 30AfiA/| - IOM2M3 + 24Aff )fc5 + 10{Mi - 3Af| - 4AfiAf3 

+12MfM2 - 6M^)k'^ + 35(A/3 - 3AfiA'/2 + 2M^)k^ + 50{M2 - Mf)k^ + 2AMik. (Al) 

I 



The ratio of the coefficients in front of the leading /c' 
terms in Kq and fq gives Hq for the Poisson distribution. 
Thus, in general, Hq differs from for a multicomponent 
Poisson distribution. 



The case of one active pair corresponds to Mr = 1, 
and the ordinary formula of NBD is restored: fq = r(A; + 
q)lT{k). 
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To demonstrate the accuracy of the wj values shown M3 = 11.0779, M4 = 34.6791, M5 = 117.238. (A2) 

in Table U for the IPPI model, we present here their more 
accurate values and moments for four active parton pairs: 

These values of Wj are larger and those of Mr are 
wi = 0.51879, W2 = 0.26914, = 0.13963, smaller than the asymptotic ones shown in Eq. 

= 0.07244, Ml = 1.76571, M2 = 4.01103, 
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